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Abstract. The X-ray Background has been resolved in the 0.5-5 keV band and found 
to consist mostly of both unabsorbed and absorbed AGN with column densities < 
lO^'^cm"^. This contrasts with the local AGN population where the column density 
range extends to Compton-thick objects and beyond (> 1.5 x lO^'^cm"^). Stacking 
analysis of the integrated emission of sources detected by XMM-Newton in the Lockman 
Hole, and by Chandra in the CDF-N and S reveals that the resolved fraction of the X- 
ray Background drops above 6 keV and is about 50% above 8 keV. The missing flux has 
the spectrum of highly absorbed AGN, making it likely that the range of column density 
at redshift one is similar to that locally, and that many AGN are as yet undetected in 
well-studied fields. 



1 Introduction 

Deep images with Chandra and XMM-Newton have shown that the X-ray Back- 
ground (XRB) is composed, at least in the 2-5 keV band, of the emission of 
many active galactic nuclei (AGN) [1] [2] [3] [4] . This was anticipated by earlier 
work with ROSAT [5] which resolved the 0.5-2 keV XRB in a similar manner. 
The XRB therefore reveals the integrated X-ray emission of AGN and is an im- 
portant tool for studying the accretion history of the Universe and the evolution 
of AGN. 

The background spectrum in the 2-10 keV band is flatter than typical AGN 
and requires that most are absorbed [6] [7] [8]. This ties in with what we know 
of AGN locally where absorbed Seyfert II galaxies outnumber Seyfert Is. What 
remains to be determined is whether powerful objects such as quasars have a 
significant obscured fraction, how the obscuration evolves with time in an object 
and how obscuration evolves with redshift within the population. 

2 Local obscured AGN 

Absorbed AGN are very common in the local Universe [9]. The three near- 
est AGN with X-ray luminosities exceeding 10"*° erg s"\ NGC4945 [10][11], the 
Circinus Galaxy [12] and Gen A [13] all lie behind absorbing column densities 
TVh > lO^^cm"^. The first two are Compton thick with A^h > 1-5 x lO^^'cm"^. 
The ratio of absorbed to unabsorbed nuclei is at least 3-5, with the column den- 
sity distribution being fairly fiat [14]. The large ratio suggests that the geometry 
is more complex than the simple unification torus model. 



2 



A.C. Fabian & M.A. Worsley 



At greater distances where higher luminosity objects are sampled there are 

some clearly highly absorbed AGN (e.g., NGC6240, [15]; IRAS09014, [16]; 3C294 
nucleus, [17]). In particular, the nuclei of many radio galaxies lie behind large 
column densities. 
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Fig. 1. The extragalactic XRB intensity as observed and the minimum absorption- 
corrected spcctnim required assuming a typical intrinsic spectrum with F = 2. 
Compton-thick obscuration or reflection could make the corrected intensity higher. 



3 Tying the evolution of AGN and the XRB in with the 
local black hole mass density 

If r] is the efficiency with which mass is turned into radiation during accretion 
(ry = 0.06 for a standard thin disc around a non-spinning black hole and a typical 
assumed value for an accreting black hole is 0.1), then 

erad(l + 2:) = mc^ 

[18] Erad is the observed energy density in that radiation now, ph is the mean mass 
density added to the black holes and z is the mean rcdshift of the population. 
Applying this formula to the intrinsic absorption-corrected spectrum of the XRB 
with a 2-10 keV bolometric correction factor of 30 and a mean source redshift 
of 2 gave a mean black hole mass density of about Ph ~ 8 x 10^ Mq Mpc'^ [19]. 

This is significantly larger than ph '~ 3 x 10^ Mq Mpc~^ found by [20] from 
the black hole mass - stellar velocity correlation applied to the local galaxy lumi- 
nosity function. Elvis et al [21] suggest that this requires the radiative efficiency 
r; > 0.15 and thus that most black holes are spinning rapidly (which boosts rj). 
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A major new result for AGN from X-ray surveys is the peak in number 

density around z ^ 0.8, not at z ^ 2 as found for quasars [22] [23] [24]. Simply 
put, the population of AGN with X-ray luminosities above 3 x 10^'* erg s~^, 
which are traditionally quasars, began dying at 2; ~ 2 whereas the population 
of less luminous objects, akin to nearby Scyfcrts, carried on to z '--^ 0.8 before 
dying. This means we should rework the Seyfert part of the Soltan equation 
with a mean z ~ 0.8. That changes the bolometric correction to 15 [17] and 
brings the predicted black hole mass density into line with observations at 4 — 
5 X lO^Mo Mpc-3 for t] ~ 0.1 [17][25][26]. 

The situation therefore looks fair for ~ 0.1 provided that all populations 
have been accounted for. A possible problem lies in the column density dis- 
tribution of local sources compared with distant ones, such as are commonly 
resolved in deep X-ray surveys. The sources found in the latter overwhelmingly 
have column densities < lO^'^ cm^^ [27] whereas locally this accounts for only 
about 40% [14]. Does the column density distribution evolve with redshift or 
have the surveys failed to find AGN with larger column densities? To study that 
we have examined just how well the spectrum of the resolved sources in deep 
X-ray fields matches the spectrum of the XRB and so whether a population of 
sources remains unresolved [28] [29] . 

4 Resolving the XRB with Deep Surveys 

The deep, pencil-beam X-ray surveys have confirmed that the XRB is made up 

of point sources. The XMM-Newton observation of the Lockman Hole (XMM- 
LH; [30]), as weh as the Chandra Deep Field North (CDF-N; [3]) and Chandra 
Deep Field South (CDF-S; [1]), are able to resolve ~ 70-90% of the XRB in the 
broad 0.5 2 keV and 2 10 keV bands. The variation in this fraction is dependent 
on the XRB normalisation^ chosen. 

These high resolved fractions are often used as the basis for claims that the 
XRB has been completely resolved up to lOkeV. Such optimistic statements 
are misleading, since it is counts in the 2-5 keV regime that dominate the 2- 
lOkeV band. In order to investigate more carefully the behaviour of the resolved 
fraction with energy we carried out a source-stacking analysis rising the XMM- 
LH data [28] . We employed a straightforward photometric approach to find the 
total resolved flux in a number of narrow energy bands. XMM-Newton was 
particularly suited to the study since it has good sensitivity up to 12keV. 

The Lockman Hole is one of the most well-studied regions of the sky and has 
received multi-wavelength attention. The unusually low level of neutral Galactic 
absorption 5 x lO^^^cm^^) makes it particularly well-suited to deep X-ray 
surveys. The total accumulated XMM-Newton exposure time has now reached ~ 
700 ks from 17 different observations. We restricted our analysis to sources within 

^ [31] summarize their own and previous work on the spectrum of the XRB which 
indicates that the normalization is probably a factor of 1.3 times greater than the 
HEA0-A2 all-sky study of [32]. See also [33]. 
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a central 10-arcmin radius with the lowest background levels and maximum 
exposure time. 

4.1 Source Detection, Photometry and Stacking 

We cleaned the raw event files to exclude background flares and bad pixels 
before extracting images in five different energy bands; 0.2-0.5 keV, 0.5-2 keV, 
2-4.5 keV, 4.5-7.5 keV and 7.5-12 keV. We carried out sliding-box source detec- 
tion to generate an initial list of bright sources. These were then masked-out 
of the images and further source-fitting carried out with maximum-likelihood 
point spread function (PSF) fitting. The final soiircc list was verified manually, 
resulting in a final fist of 126 detected sources located within 10' of the observa- 
tion. (Since this early work our subsequent re-analysis has made use of a more 
complete source list of 156 objects although the additional 30 sources make no 
significant difference to our results.) 

We adopted straightforward aperture photometry to determine the count- 
rate of each source and in each energy band. The three XMM-Newton instru- 
ments were processed separately to provide robustness and because of differences 
in background level and PSF. For each instrument and energy band an exposure- 
corrected image was made by dividing by the exposure map. In each case we also 
created a source-free image by masking out all the detected sources. For each 
source and in each energy band, the total count-rate was extracted from within a 
circular aperture centred on the source position. The local background level was 
measured in an annulus surrounding the source but here we used the source-free 
version of the images to avoid contamination. Given the aperture-measured and 
anmilus-measured count-rates the true source and background count-rates can 
be easily determined. 

When calculating the true source and background count-rates we assumed 
an analytical fit to the PSF. Tests of the quality of the analytical model on a 
number of sources found it to be acceptable within the energy bands and off- 
axis angles required. The radius of the circular extraction aperture was chosen 
to maximise signal-to-noise in the final source count-rate. The background- and 
exposure-corrected count-rates for each source were converted to fluxes using 
conversion factors computed with XSPEC assuming a F = 1.4 source plus A^h = 
5 X 10^^ cm~^ Galactic absorption. 

The total resolved flux in each energy band was calculated by summing the 
measured flux for each source. Every source is included regardless of whether 
or not it is detected in the band. Flux errors were added appropriately and 
combined with the estimated error in the counts-to-flux conversion. A correction 
for Galactic absorption was made; this is ~ 15% in the 0.2-0.5 keV band and 
< 1% in the harder bands. 

4.2 XMM-Newton Results 

Figure 2 shows the resolved XRB intensity for each of the XMM-Newton instru- 
ments along with the total XRB. Prom 0.5-12 keV we took the total XRB to be a 
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r = 1.41 power-law with a normalisation of 11.6 keV cm~^ sr~^ keV~^ [31]. 
In the softest band, 0.2-0.5 keV the extragalactic background level is uncertain 
and is estimated to lie at 20-35 keV cm~^ s~^ sr~^ keV~^ [34]. 
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Fig. 2. The extragalactic XRB intensity resolved from detected sources in the energy 
bands. The values from the three instruments are plotted, for clarity, as points in the 
centre of each band and offset horizontally with respect to each other. Squares, circles 
and triangles represent the PN, MOS-1 and MOS-2 instruments respectively. Errors 
are one sigma. The bars represent the total extragalactic XRB intensity. 



From 0.2-2 keV the resolved fraction of the background is consistent with 
~ 90%. The 2-4.5 keV band also shows a high resolved fraction of 70-80% but the 
harder bands > 4.5 keV only resolve < 60%. The missing XRB component has a 
spectral signature consistent with what would be expected from an undetected 
population of highly obscured sources. 

The level of the unresolved fraction depends upon the Lockman Hole being 
a representative sample of the sky, as well as the normalisation taken for the 
total XRB, which itself is still subject to errors of 10-15%. We do not correct 
for the bright-end source population: the pencil-beam nature of the field means 
that bright sources 5 x 10""'^^ erg cm~^ s^^) are not well-sampled whilst those 
brighter than ~ 10~^^ erg cm~^ are not seen at all. The additional contribu- 
tion by these bright objects would increase the resolved fraction by 10-20%, 
but this is only comparable to the effect of field-to-field variations. Critically, 
since bright sources typically have very soft spectral slopes, no additional contri- 
bution can account for the missing XRB above 5 keV, nor could any adjustments 
to XRB normalisation. 
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4.3 The Chandra Deep Fields 

Although Chandra has low effective area at energies exceeding 7 keV, the CDF- 
N and CDF-S probe much deeper than the XMM-LH. Wc used the CDFs to 
assess whether the contribution from the fainter sources detected by Chandra 
were able to explain the XMM-Newton results. To aid comparison, as well as to 
take advantage of numerous software improvements and the complete XMM-LH 
source catalogue, we re-analysed the XMM-Newton data using identical energy 
bands to the Chandra work (see [29]). 

The CDF-N and CDF-S cover 447.8 and 391.3 arcmin^ respectively Wc used 
the catalogues from [3] which contain 503 and 326 sources respectively for the 
fields as well as photometry. Source-extraction apertures enclosing 90-100% of 
the counts were used and the local background correction was determined from 
Poisson statistics. Counts-to-flux conversion was performed, where possible, us- 
ing conversion factors appropriate to the hardness ratio of the source. For sources 
with an insuflacient number of counts to permit this, a T = 1.4 spectrum was 
assumed. 

Prior to our stacking analysis we applied corrections for Galactic absorption 
as well as several systematic effects which have been recently quantified [35]. 
The total resolved flux in each of five photometric bands (0.5-1 keV, 1-2 keV, 
2-4 keV and 4-6 keV) was calculated by summing the measured flux for each 
detected source. In order to calculate a resolved intensity the flux from each 
source must be divided by a solid angle on the sky; however, since the CDFs 
show a substantial increase in sensitivity towards the centre of the fields, this 
solid angle is a function of source brightness. The total solid angles of 447.8 
and 391.3 arcmin^ are only applicable to the brightest objects, whereas fainter 
sources are only detectable over a fraction of these regions and the appropriate 
solid angle to take is smaller. We use the sensitivity functions from [3] when 
calculating source intensities. 

4.4 Bright-End Correction 

In order to remove the 'noise' in the resolved fluxes that is due to poor sampling 
of the brighter sources (which only occur with sky densities of ~ a few on the 
size of the deep surveys), we ignored sources with 0.5-8 keV fluxes greater than 
5 X 10~^^erg cm~^ s~^. We then used number counts from wide-area surveys 

to replace the missing flux due to all the sources brighter than this. Wc used 
the the logA^-logS' functions from [36], compiled in the 0.5-2 keV and 2-10 keV 
bands. 

A number of studies have confirmed a strong dependence of source spectral 
shape with flux. At the bright-end sources have spectral slopes much softer than 
the r' = 1.4 of the XRB. We took a linear flt to the relationship as observed by 
[37], although we conservatively impose a maximum of r' = 2. Given this relation 
we can compute the bright-source contribution for each energy band, using the 
appropriate spectral index at each source flux and integrating the logA^'-logS' 
function. 
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4.5 Results 

Figure 3 shows the bright-end corrected and resolved XRB as seen in the CDF-N, 
CDF-S and XMM-LH. The extragalactic XRB level is taken to be the spectrum 
measured by [31], but at high energies we take account of the slight turn-over, 
and below 1 keV we take account of the steeping over the spectrum implied by 
AGN-shadowing observations (for details refer to [29]). 

The resolved fraction is ~ 70-90% up to ~ 4keV. XMM-Newton sees a 
decrease in the 4-6 keV band whereas the CDFs remain consistent with a high 
resolved fraction of ^ 80-90%. This difference is due to the faint, hard sources 
which are detected in the CDFs but not the (shallower) XMM-LH. In the 6- 
8keV range all three surveys only resolve 50-70%, whilst in the 8-12 keV band 
XMM-Newton only resolves 40-60% of the background. The total resolved flux 
in the broad 0.5-8 range remains high (80-85%), and illustrates the need for 
narrow-band stacking to reveal the decrease > 6 keV. 
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Fig. 3. The resolved fraction of the XRB as seen in the CDF-N (black crosses), CDF- 
S (grey crosses) and XMM-LH (the combined PN/MOS-l/MOS-2 result; black dia^ 
monds), as a function of energy. 



4.6 The Missing Source Population 

We have modelled the spectral shape of the missing XRB fraction as a single 
population of sources. We assumed an underlying F = 2 power-law, inclusive of 
an i? = 1 reflection component, plus intrinsic photoelectric absorption. A grid 
of spectra was constructed over the range z = 0.1-3 and iVn = 10^^-10^^ cm~^ 
and for each we calculated the goodness-of-fit between the model spectrum and 
the unresolved XRB spectrum. The confidence contours in source spectral shape 
are shown in Figure 4. 
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Fig. 4. Confidence contours (horizontal) of fits between the missing XRB spectrum 
and an obscured power-law model for the XMM-Ne*vvton data. 68, 90 and 95% contours 
are shown. The best-fitting model has z ~ 0.8 and ATh ~ 4.5 x 10^^ cm~^. The other 
contours indicate the maximum, unobscured, rest-frame 2-10 keV luminosity of a single 
source in the population. 

The XMM-LH and CDFs show similar contours for the characteristics of 
the missing population, although it is the XMM-Newton results, -which extend 
up to 12keV, which provide the best constraints. The best-fitting models have 
z = 0.5 1.5 and A^h ~ 10^^-10^'* cm^^, pointing to a highly obscured population 
of AGN. Since the missing sources are undetected, the survey sensitivity allows 
an upper limit to be placed on the luminosity of each object. For the best-fitting 
source models this upper limit is < 5 x 10''^ erg s~^. The sources must occur with 
areal densities of > 2800 deg~^ - a considerable number. This would correspond 
to at least ~ 350 undetected AGN in the CDF-N, and suggests that there must 
be > 3 times more undetected, obscured sources than detected, unobscured ones. 
Figure 5 shows the resolved and total XRB with the typical shape of the sources 
required to account for the missing intensity. 

The absorption-corrected flux of the missing population can be used to es- 
timate the black hole mass density of the obscured population. Assuming an 
accretion efficiency of 0.1, and a bolometric correction of 0.1, the missing popu- 
lation would correspond to a black hole mass density of ~ 0.5-1 x 10^ Mpc"'^; 
~ 10-20% of the total local value. 

5 Summary 

Studies of the XRB continue to reveal the evolution and properties of AGN. 
Most of the XRB is resolved up to an energy of 5 keV and is due to AGN - 
Seyferts and quasars - with absorbing column densities < 2 x lO^^cm"^. The 



The Obscured X-ray Background and Evolution of AGN 



9 




0.1 1 10 100 

Energy (keV) 



Fig. 5. The total extragalactic XRB (upper solid curve) along with that resolved in 
the XMM-LH (black diamonds), CDF-N (black crosses) and CDF-S (grey crosses). The 
lower solid curve shows the spectrum of the best-fitting model for the missing source 
population at 2 ~ 0.8 and with Nu ~ 4.5 x 10^^ cm~^. 

identification of large samples of such objects shows that the lower luminosity 

ones (Seyfert-like) evolve rapidly to z ^ 0.8 while quasars evolve more slowly to 
z ~ 2. Assuming that the Seyfert-like ones have lower black hole masses, those 
black holes continued to grow at later stages of the Universe than the more 
massive ones associated with quasars. 

The resolved fraction of the XRB drops to 60-70% above 6 keV and 50% 
above 8 keV. The missing fraction between 6-8 keV is not due to the known 
faint sources in the 2-5 kcV band since it would have been detected by the 
stacking analysis which corresponds to an exposure time of 0.5 Gs or more. The 
missing fraction is most simply explained as due to the high column density 
(A^H > 2 X lO^'^cm"^) extension to the observed Scyfcrt-likc objects evolving 
rapidly to z 1. Such a population was implicit in XRB synthesis models such 
as that of [38] or [39]. 

Presumably there is also a Compton-thick population as well, which would 
make less of an impact on the XRB since Compton down-scattering then depletes 
the emergent X-ray luminosity. The total black hole density due to Compton- 
thin sources is then about w 5 x lO^rjZl Mq Mpc^^, where rj = 0.1?7_i. This 
is consistent with recent estimates of the local value of ph [25] [26] for rj = 0.1. 

The issue of highly obscured and Compton-thick quasars remains unresolved 
[40]. Powerful radio galaxies represent one, relatively small, population. Some 
powerful, highly obscured, radio-quiet quasars have been found [41] [42] [43] [44], 
with [45] estimating a surface density of at least 10 deg~^, The uncertainties are 
currently too large to assess their general importance. Future work combining 
X-ray and Spitzer data will be instructive. 
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